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SECTION  I 
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I 

INTRODUCTION 


In  a companion  report  [l],  the  static  analysis  of  a 
conical  antenna  over  a ground  plane  is  presented.  In  this 
report,  the  analysis  is  extended  to  treat  the  time-har- 
monic case  and  to  incorporate  a model  of  the  resistive 
loading  of  the  structure.  The  resistive  loading  is  intended 
to  reduce  the  effect  of  diffraction  from  the  cone  edge  at 
the  higher  frequencies. 

If  the  cone  has  no  topcap,  the  analysis  may  be  consi- 
derably simplified  and  a simple  integral  equation  for  this 
situation  is  derived  in  Section  II  and  implementation  of 
a moment  method  solution  is  considered  in  Section  III.  In 
Section  IV,  an  integral  equation  for  the  cone  with  a top- 
cap  is  derived.  Presented  in  Section  V are  numerical  results 
in  the  frequency  domain  for  currents  on  a loaded  conical  antenna 
both  with  and  without  a topcap.  In  Appendix  A expressions 
are  derived  for  the  computation  of  fields  from  the  currents 
and  Appendix  B gives  the  derivation  of  an  alternate  integral 
equation  from  that  derived  in  Section  IV. 
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SECTION  II 


FORMULATION  OF  AN  INTEGRAL  EQUATION  FOR 
A BICONE  WITHOUT  ENDCAPS 

For  the  symmetrically  driven  biconical  structure  of 

Figure  1,  the  current  on  both  the  cone  and  its  image  are 

radially  directed  and  have  no  circumferential  (0)  variation. 

Hence  the  magnetic  field  tangent  to  the  cone  is  ij)-directed 

and  the  boundary  conditions  can  be  satisfied  by  fields  which 

are  '■ransverse  magnetic  (TM)  to  r.  Thus,  the  fields  may 

npletely  determined  by  a radially-directed  vector 

ncial  A * A^r  [2] . In  an  eigenfunction  solution  to 

such  problems,  the  fields  are  determined  in  the  bicone 

region  from  a vector  potential  A^  which  comes  from  a 

homogeneous  solution  of  the  wave  equation.  In  order  to 

derive  an  integral  equation,  however,  A^  must  be  expressed 

in  terms  of  the  current  on  the  bicone.  In  particular,  a 

free  space  Green's  function  is  to  be  found  for  the  vector 

potential  A^  due  to  a unit  radially-directed  current  element. 

A superposition  integral  then  expresses  the  total  vector 

potential  due  to  currents  on  the  cone. 

Beginning  with  the  assumption  that  the  magnetic  field 

is  determined  from  A = A f, 

r 

H = - VxA 
U 

2 

I 


Figure  1.  Ceonietry  of  cone  over  a ground  plane. 
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and  using  Maxwell's  equations. 


VxE  = -jujyH 
VxH  = jo)eE+J 


one  readily  determines  by  standard  procedures  that  A satis- 
fies the  vector  Helmholtz  equation 
- - 

VxVxA  - k“A  = uJ^r  - ja)UeV$  (1) 

where  <1>  is  a scalar  such  that 

E = -jo)A  - V<1> 


For  a radially-directed  unit  current  element, 

j = J e = 

r 

- <5(r-r')6(6-e')6((j)^’) 

r ' ^ sin  9 ' 


Expanding  out  (1)  yields 


-1 


2 . a 39  i 

r s 1 n 9 ^ 


[ 3A 


2 , 2- 
r sin  9 


3^A 


9 

k“A 


1 - 

r 3r99j  ^ [ r sin  9 9<])3r  J 


9 

1 ^“^r 


pj^r  - jojye 


9<J 


1 94 
7 9^ 


1 


r sin 


94  : 

~ *TT“  ^ 


(3) 
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from  which  it  is  seen  that 


3^A 
r 

3 r 3 0 


3‘I> 

■jujue  ^ 


3^A 
r 

3 4i3  r 


3'P 


(4) 


These  conditions  are  automatically  satisfied  by  the  gauge 
choice 


$ 


1 


3 A 


jujue 


3r 


(5) 


Substituting  (2)  and  (5)  into  (3)  leaves  t'*'-  >ne  scalar 
component  equality 


3~A 


3 r 


2 . „ 39 

r s in  9 


sin  o 


3A 

I 

30 


2 . 2, 

r s in  d 


9 

3“A 


3 (S  ‘ 


+ k“A  = -u6(r-r')  (6) 

r 


which  can  be  rewritten  in  the  more  convenient  form. 


(V^+  k^)(  — ) 


-u6  (r-r  ' ) 
r ' 


To  obtain  (7),  one  notes  that 


5(r-r  ' ) 
r 


6 ( r - r ' ) 
r ' 


(7) 


A solution  of  (7)  which  satisfies  the  radiation  condition  for 
a exp(j'Jt)  time  convention  may  be  written  by  inspection  of  (7) 
as 


5 


A = 


U r e 


-jklF-;' 


r Attt  ' r-f ' 


(8) 


and  the  general  solution  to  (3)  for  a distributed  set  of 
currents  is  obtained  by  superposition: 


r Att 


I -jk|F-r! 

I I 

/ r'lr-r'l 

•/ V 


dV  ' 


(9) 


This  forni  of  the  vector  potential  has  also  been  used  by 
others  [3].  For  the  syinme  t r i ca  1 ly -ex  c 1 1 ed  cone  and  its 
image , 

•J„  ^ ) r 1 

(10) 


J^(r' ) 


J ( r ’ ) p 

' <S(0'-0g)  - 6(0'-TT+eQ)J 


' *•  V r 
L£ 5 

r ' 


where  J is  the  bicone  surface  current  density.  Sub- 
s r 

stituting  (10)  into  (9)  gives  the  desired  equation  for  the 
vector  potential: 


2tt  L 


Ur  sin  0 ^ 
A TT 


// 

0 0 


J (r  ' ) 

s r 


^-jkR  ^-jkR 


dr ' d0 ’ 

(11) 


whe  re 


r“  + r 


' 2 _ 


rr'  [sin  0 sin  0 -cos  (;j)-0  ' )+  cos  0 cos  0 


0- 


The  plus  superscript  denotes  source  points  on  the  upper  bicone 
surface  while  the  negative  sign  denotes  source  points  on  the 
image  surface. 


6 


,n 

I 


It  is  convenient  to  introduce  the  total  axial  current 


I(r')  » 2nr’  sin  0_J  (r') 

0 sr 


(12) 


so  that  (11)  becomes 

2TT  L 

/ 

A = 


ur  r r Kr' 

8tt-  J J r’ 


g-jkR  ^-jkR 


dr'd$'  (13) 


0 0 


The  radial  component  of  electric  field  is  now  given  by 


E 

r 


L_ 

j o)ue 


( 14) 


The  simplicity  of  (13)  and  (14)  compared  to  the  usual 
vector  potential  representations  should  be  emphasized  at 
this  point.  One  notes  that  in  the  usual  representation,  two 
vector  potential  components,  and  Ag , would  be  present. 
Furthermore,  the  integrands  of  the  potential  integrals  would 
contain  somewhat  complicated  dependences  on  angles  between 
observation  and  source  points  which  arise  from  projecting 
the  source  vector  onto  the  potential  component  vector  for 
each  source  and  observation  point.  Finally,  the  expression 
for  the  radially-directed  electric  field  would  be  complicated 
and  difficult  to  handle  numerically  compared  to  the  approach 
to  be  followed  here.  These  complications  indeed  will  appear 
in  the  formulation  which  includes  a topcap  on  the  bicone  struc- 
ture (Appendix  B). 
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An  integral  equation  for  the  current  is  obtained  by 
applying  the  boundary  condition  that  the  radial  electric 
field  must  equal  the  impedance  loading  times  the  total 
current  density,  i.e., 

> Z^(r)I(r) 

Since  all  currents  and  fields  are  ^-independent  , it  suffices 
to  take  all  observation  points  along  the  intersection  of  the 
plane  = 0 and  the  conical  surface.  Hence,  we  obtain  finally 


1 

ju)ue 


A - 2 (r) I(r) 
r s 


= 0,  0<T^l,  9 = 9q,'$  = 0 

(15) 


Equation  (15)  is  an  in t eg r o -d i f f e r en t ial  equation  for  the 
induced  current  on  the  bicone.  As  it  stands,  (15)  does  not 
appear  to  contain  a driving  term  due  to  the  applied 
voltage  at  the  bicone  terminals.  In  the  next  section, 
however,  this  term  appears  as  a "boundary"  condition  on 
dA^/dr  at  r=0.  One  also  notes  in  (15)  that  discrete  or 
lumped  loading  may  be  introduced  by  allowing  Z^(r)  to  be 
represented  by  appropriate  6-functions, 

n*  1 

for  loads  where  Z^^^  is  the  impedance  of  the  n^^  load 

loca  ted  at  r = r . 

Ln 
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SECTION  III 

APPLICATION  OF  METHOD  OF  MOMENTS  ID  A 
BICONE  WITHOUT  ENDCAPS 

The  usual  procedure  in  applying  the  method  of  moments 
[4]  is  to  first  represent  the  unknown  current  as  a linear 
combination  of  an  appropriate  set  of  basis  functions  and 
then  "test"  the  resulting  integral  equation  with  a series 
of  testing  functions.  Here  it  is  convenient  to  reverse 
this  order  and  to  first  test  the  equation  before  expanding 
the  current.  A set  of  testing  functions  which  offer  a 
number  of  advantages  in  a numerical  procedure  are  the  piece- 
wise  sinusoidal  testing  functions: 


sin  k(Ar-r) 
sin  kAr 

0 


0 ^ r A r 
Ar  , < r < L 


where  Ar  » L/M,  r » (m-l)Ar,  m«l,2,...M 

m 

These  testing  functions  are  shown  in  Figure  2.  .An  inner 
product  is  next  detined  as 


9 


< f ( r ) , g ( r ) > 


(17) 


/ 


f ( r ) g (r ) d r 


and  (15)  is  successively  tested  with  each  of  the  w^, 
ni»l,2 M: 


dr- 


A , w 
r m. 


Z (r) I(r)  , w 
s m , 


0,  m-1,2, . . . ,M 


(18) 


Taking,  for  ihe  moment,  m=l  and  integrating  the  first  term 
by  parts  twice  results  in 


A , 

dA 

r 

i * - — — 

I 2 

! dr 

r 

dr 

r = 0 


sin  kAr  r'"2 


A, (r_  ) 


k cos  kAr 


A_  (r»  ) 


whe  re  A(r)=*A 
r m r 


r-r 


9 = 9, 


sin  kAr  "r"'2' 
. Note  that  although 


9-TT/2 


A^(tj^)  is  not  zero  along  the  bicone 


, 3^A 

1 r 


0. 


(19) 


In  fact,  one  notes  that 


9 jojper  3r39 

and  that  the  bicone  voltage  at  r=0  is  just 


11 


tt/  2 


/ 


’'o-  ; 


r="0 


>1/2 


j wpe  I 


3^A 
r 

3r39 


d0 


1 


r = 0 


j u)U  e 3 r 


where,  using  (19),  one  sees  that  3A^/3r*0  at  Q^tt/; 
Thus  for  tn  = 1 , (18)  becomes 


jojvie  sin  kAr 


[- 


cos  kAr  A^(rj^)  + A^(r2 


>] 


Zg(r)I(r)  ,wi^=  -Vq 


(20) 


For  m=  2,3,i,...,M,  integration  by  parts  twice  in  (18)  re- 
sults in 


^ r-7 Fa  (r  ) 

jojue  sin  kAr  r mfl 

m « 2 , 3 , . . . ,M 

(21) 

Note  that  the  choice  of  testing  functions  has  resulted 
in  removing  all  the  derivative  operations  from  the  operator 
equations.  This  is  the  principal  advantage  of  the  testing 
functions  chosen. 


2 cos  kArA  (r  ) + A 
r m 


r^"m-l^] 


Z (r)I(r),  w 0 

® m . 
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A matrix  equation  now  results  if  the  current  is 


expanded  in  an  appropriate  set  of  basis  functions.  A 
convenient  set  is  the  pulse  functions  defined  by 


(r) 


r 1 , 0 < r < Ar/2 

lo  , Ar/2^r^L 


Pn(r) 


r 1 . I r-r^  I <_  Ar/2 
0 , I r-r^l  > Ar/2 

n - 2 . 3 , . . , M 


(see  Figure  3)  and  the  resulting  current  expansion  is 


I (r  ) 


M 


r 

= I 

n»  1 


I_P„(r) 
n n 


(22) 


Note  that  the  current  at  the  bicone  edge  r«L  is  auto- 
matically zero  by  our  choice  of  basis  functions  (Figure  3). 

When  (22)  is  substituted  into  (20)  and  (21),  there 
results  the  system  of  linear  equations 


jwuG  sin  kAr 


[- 


cos  kAr  l'(rj^;rj^,rj  + ) * f(r 


2 ’ ^ 2 ’ 1 U 


< 


Z(r)p(r).  w 
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M 

^ I I 

n»  2 


n \ jdJUG  sin  k^r 


cos  kAr  y (r  ; r^_  , r^+)  + y (r  ^ ; r , r^+)J 


Z (r)p^(r),  w >)=  -V 


s n 


0 

(23) 


and 


lAjoJUe  sin  kAr 


+ y(r 


m+1 


;ri.ri+^ 


<\^(r)Pi(r).  ^ 


^ n I 1 (jju e sin  kAr 
n * 2 I 


\r['^^"m-l*"n-’"n+^-  ^ "m  ’ "n“  ’ "n+^ 


+ y(r 


m+1 


’"n-*^  + 3 


/ 


Z„(r)p_(r),  w 


s n 


0 ,m-2 , 3 , . . ,M 

(24) 


These  equations  may  be  assembled  into  the  matrix 


equat ion 


whe  r e 


ZI  - V 


“ — 

”-v 

I 

0 

71 

0 

* 

, V - 

0 

Iw 

• 

M 

0 
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and  the  elements  of  the  impedance  matrix  may  identified 

from  (23)  and  (24).  The  functions  4' ( r ; r _,r  +)  are 

n n 

def ined  by 


(25) 


where 


R » r^+  r'^  - 2rr'(sin^0Q  cos  ())'  +cos^0q)  (26) 

and  where 


r ± = r + Ar/2 
n n — 

The  evaluation  of  the  double  integral  in  ( 24)  is 
simplified  by  analytically  approximating  the  integration 
with  respect  to  r'.  This  may  be  accomplished  by  noting  that 
in  general  k(r^+-  r^_)  <<  1,  so  that  a few  terms  in  a 
Taylor  series  expansion  about  some  point  r^  in  the  interval 
fr  ,r  j.]  should  be  accurate.  Accord  ing  Iv,  one  writes 


-jkR-  'Jk(R--R^) 

e ” e e 

= e " [1-j  k(R“-  R-) ] (27) 


where  R « R . Substituting  (26)  into  (24)  and  listing 

n I r * r 

n 

as  a fourth  argument  the  point  about  which  the  expansion  is 


made,  one  obtains 
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^ r ; r ) 

2^  , 


47T 


I^J  /■[• 


— jkR"*"  fl+jlcR"*”  ., 
n ! -*  n .1  k 


i r ' r'*'  r'^ 


0 r 


-jkR, 


r ' R 
1+jkR 


r ' R 


s_  _ ^ ! 


r j 


d r ' 


’T  f , 

1 

r 

1+jkR 

1 • 

n , 

J.  log 

0 1 

- 

r +(R  + r - r b ) 

n n- n- 

r (R"*”.  + r - r .b"'’) 
n-  n+  n+ 


■j  kiln 


r + 
n 


-jkR, 


1+jkR 


log 


r . (R  _ + r 
n 


r_(R.+r-r,b) 
n n+  n + 


- j k£n 


n+ 


r _ 
n _J 


d<t>  ' 


(28) 


±2  2 

where  b » sin  0^  cos  4> ' + cos  9^.  In  several  situations, 
appropriate  limits  of  the  integrand  of  (27)  need  to  be  taken 
First,  when  the  source  is  the  current  segment  at  the  bicone 
terminals,  the  integral  (27)  reduces  to  the  simple  form 


1'  ( r ; r ^ , r j^  + ) 


'f  (r  ; r^  , r^+,  r ^) 


II 


-j  kr 

e ( 1+j  kr ) log 


^1+  + 


b r + 


R]_+  + r 


b''’r 


9 + 


dil ' 


(29) 
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As  the  observation  point  r in  (28)  approaches  the 


bicone  terminals,  r-*-r^=  0,  the  limiting  form  of  the  inte- 
grand can  be  integrated.  The  result  is 


nri;ri, 


1^1+, ri> 


u 

21T 


log (cot 


(30) 


a very  interesting  result  that  is  independent  of  the  sub- 

domain  size  at  the  bicone  terminals.  Finally,  all  the  so- 

called  "self  terms"  4'(r  ;r  ,r  , ,r  ) . n^^l,  contain  an 

n n-  ’ n+  n 

integrable  singularity.  In  fact  one  easily  establishes 
that 


log 


r _^(r"^  + r - b'^r  ) 
n-b  n- n_2 

r .(R"*”,  + r - b'*'r  ,) 

n n+  n+ 


r-<-r 


- 2 log  1 $ ' 


This  singular  term  is  then  subtracted  from  the  integrand  in 
(27),  resulting  in  a non-singular  integrand  which  is  then 
numerically  integrated,  The  term 


TT 


/ 


Jin  I 0 ' I d p ' 


(JlnTT  -1) 


is  then  added  to  the  result  to  take  care  of  the  part  of  the 
integral  contributed  by  the  singularity. 
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SECTION  IV 


FORMULATION  AND  NUMERICAL  SOLUTION  OF  AN 
INTEGRAL  EQUATION  FOR  A BICONE  WITH  ENDCAPS 

The  formulation  of  the  integral  equation  for  a cone 
radiator  over  a ground  plane  with  an  endcap  is  considerably 
more  complicated  than  that  for  the  case  when  the  endcap  is 
not  present.  It  is  possible,  however,  to  generalize  the 
approach  used  for  the  bicone  without  a topcap  and  to  trans- 
form the  derivatives  appearing  in  the  equations  into  har- 
monic operators  along  the  radial  cone  and  topcap  coordinates, 
as  is  done  in  Appendix  B . This  approach  has  the  advantage 
again  that  testing  with  piecewise  sinusoids  allows  the  re- 
placement of  derivatives  by  a finite  difference  of  potentials. 
However,  to  effect  this  transformation,  an  extremely  com- 
plicated kernel  must  be  used  (see  Appendix  B ) which  contains 
many  singularities  other  than  the  usual  ones  where  source  and 
field  points  coincide.  While  this  approach  has  been  used, 
it  has  been  found  to  be  unwieldy  and  rather  inefficient. 

The  approach  described  here  begins  with  the  description 
of  fields  in  terms  of  the  more  commonly  used  vector  magnetic 
and  scalar  potentials  expressed  in  terms  of  the  cone  currents 
and  charge.  However,  it  is  found  that  these  potentials  are 
singular  at  the  bicone  terminals  which  again  creates  an 
unnecessary  complication.  In  order  to  circumvent  this  problem, 
the  cone  and  Image  surfaces  are  allowed  to  intersect  with  a 
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small  "waist"  of  radius  "a"  (Figure  4).  If  "a"  is  very 
small,  there  should  be  negligible  difference  in  the  input 
impedance  and  currents  found  for  this  case  and  that  for  the 
limiting  case  of  a=0.  For  convenience,  the  cone  coordinates 
are  defined  with  respect  to  the  projection  of  the  cone  sur- 
face to  a tip,  as  in  Figure  4.  Furthermore,  the  direction 
of  the  unit  vector  r^  and  the  positive  direction  of  corres- 
ponding vector  components  is  taken  to  be  towards  the  center 
of  the  topcap,  in  the  direction  of  decreasing  r^. 

The  integral  equations  are  obtained  by  setting  the 
radiated  field  tangent  to  the  cone  surface  equal  to  the 
impedance  drop  per  unit  length  due  to  the  loading: 

-jcjA  - 1-^  - Z I =0,  a/sin  r < L + a/sin  0_  (31) 

■^rdr  sc  Oc—  0 

c c 


3 <5 

-jwA  - T — - Z = 0,  0<  r < L sin  9_  + a (32) 

r dr^  st  t u 

t t 

where  the  tangential  components  of  magnetic  vector  potential, 

A and  A , are  given  by 
"^c  t 


■^p 


2tt  L+a  /sin  0 
I 


II 


D -jkR  -jkR 

0|-  j.+  '^pc  j.--'pc 

coste  cos^e 




d r ' d({) 
c 


0 a/ sin  0, 


2tt  Lsin0_  + a 

0 r 


II 


cos 

LL 


+ -jkR 


pc 

+ 

Pt 


pc 


-jkR 


cos  C ^ e 
+ ^ 


pt- 


dr  ^d-?  ' 


pt 


pt 


c , t 


(33> 
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The  scalar  potential  is  given  by 


-1 


8tt~  j u)E 


2tt  L+a / sin  9 


/ I 


dl 


c / e 


-ikR  -jkR 

pc  „ 


pc 


d r 


dr  d(J)  ' 
c 


0 a / s in  9 , 


pc 


pc 


2tt  k sin  0Q+a 


/ / 


dl 


t e 


-jkRpt  -jkR 


pC 


d r 


dr  d4>  ' 
t 


Pt 


Pt 


p = c , t (34) 


The  currents  I and  I are  the  total  linear  currents  on 
c t 

the  conical  and  topcap  surfaces,  respectively,  and  are  re- 
lated to  the  corresponding  surface  current  densities  and 


J by 
t 


I = 27rr  sin  9,.  J 
c c 0 c 


I = 2irr  J 
t t t 


(35) 


The  distance  quantities  are  all  of  the  form 


* / ' 2 ± ' ± 

R » r +2br+c  , p,q  = c,t 

pq  q Pq  q pq 


where 

* 2 2 - 
b “ « - r sin  9„  cos  0'  + r cos  9 - a cot  9-cos  9^(1+1) 

cccOCq  00 


c ■ * + a^coC  2r  a cos  9„cot  9^(l-^i) 

cc  c 0 c 0 0 

b ~ =-r  sin  9„  cos  li ' 
ct  c 0 

t 2 1 

c = r - 2r  cos  6»(a  cot  9_+  L cos  9_)  + (a  cot  0„+  L cos  9„) 
etc  c U 0—  0 0—  0 

+ - 7 

b = + L cos"9-  - a cos  9_cot  6„  - sin  9~  cos  6' 
tc  0 0 0 t 0 

±2  2 
c^^  = r^+(L  cos  00+3  cot  9q) 

+ 

b^^  = -r  cos  d)' 

1 1 c ^ 

+ ? 2 2 ~ 2 
+ L cos  0q(1  + 1) 

The  angles  between  the  source  current  elements  and  the  tan- 
gential component  of  electric  field  at  the  observation  point 
are  determined  by 

- 2 2 
cos  5 » + cos  (J)  ' sin  0^  + cos  0. 

^cc  _ 0 0 

+ _ 

COS  5 * + sin  0^  cos  $ * 

c c 0 

+ 

cos  * i cos  il  ’ 

cos  5 ” “ + cos  <p  ' sin  0- 
t c u 
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It  Is  convenient  to  choose  as  testing  functions  the 
pulse  functions  shown  in  Figure  5.  Thus,  testing  (31)  with 
p^  results  in 


■l’l>  ■ 0 


Upon  integrating  by  parts  in  the  central  term,  and  noting 

that  A is  slowly  varying  over  the  interval  and  hence  may 
c 

be  approximated  by  obtains 

c 


Ar 

-ji^Ar  (r^^) 
c 


- Ar^/2)  - 


But  Is  just  the  bicone  terminal  voltage  Vq  with  res- 

pect to  the  ground  plane.  Hence, 


-j'^Ar  (r^^)Ar^  - 2<t>(r^^+  Ar^/2)-  2<^Z^I^  , p - -2V^  (36) 

For  the  remaining  testing  functions  on  the  cone, 
testing  of  (31),  Integration  by  parts  on  the  scalar  potential 
term  and  approximation  of  the  vector  potential  by  its  value 
at  the  center  of  the  pulse  yields 


j'ijArA  (r  )-['t'(r  +Ar/2)-5(r  -Ar/o) 

•’  c r Cm  ‘ cm  c cm  c' 

c 
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(37) 


• 1 


At  the  edge,  the  testing  pulse  consists  of  two  parts, 
one  on  the  cone  surface  and  one  on  the  topcap.  Hence,  both 
equations  (31''  and  (32)  must  be  used.  With  approximations  on 
the  two  vector  potential  components  similar  to  that  above, 
integration  by  parts,  and  enforcement  of  continuity  of 
the  scalar  potential  at  the  edge,  one  obtains 

rAr  Ar  'v.  -1 

-jo)  (r^,^  ) + -5-  A 

c c t 

-Ar^/2)-$(r^^+Ar^/2)]  - ^ = 0 (38) 

‘ c ' c ' 

where  I = I or  I as  is  appropriate.  On  the  topcap, 

r c t 

one  has,  analogous  to  (37), 


(r^^)  - ic^/2)  - Ar^/2)1 


- <^,’t  • P»  ■ »•  "-2 <”> 


The  current  is  next  expanded  in  the  set  of  pulse  functions 


p of  Figure  5, 
n 


N 

c 

I (r  ) = y I p„(r  ) 
c c nn  c 

n- 1 

't  (■^t>  ’ ^Pn^'t^ 

n*N 

c 


(^0) 


(^  1) 
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Note  chat  the  current  at  the  edge  of  Che  cone  is  the 

c 

same  on  both  the  cone  and  the  topcap  surfaces.  The  derivatives 
of  the  currents  above  are  approximated  by  a finite  difference 
of  adjacent  current  pulses  which  is  then  assumed  to  be  ex- 
panded in  its  own  set  of  "charge"  pulses  (see  Figure  6 ); 


dl 

C r 

dr  ' ^ 
c n«  1 


"n+1  ^n 


Ar 


ax  N fN 
dl  ^ c t I 


d r 


- = I 

t n = N 


n + 1 n 


Ar , 


(42) 


(43) 


where  I is  taken  to  be  zero. 

^c  ' t 

Thus  the  vector  potential  quantities  in  (36)-(39>  may 
be  written  as 


cl 


r , + Ar  /2 ) 
cl  c 


r + Ar  /2 ) 
cn  c 


Iv  ['i'  (r  .r  „ - Ar  /2,  r „ ) + 'f  (r  ,r  .,r  .+  r /2)] 
*■  pc  p cN^  c cN^  pt  p tl  tl  t 


N +N^ 
c t 


1 I 'f  (r  , r 

^ n n h n 


n«N  +1 
c 


n p t p t 


.ntl-N  - ^"t^2  ’ "t,n  + l-N  ^ ^^/^)  j 
c ll  J 


p » c , t 


(44) 
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and  the  scalar  potential  is 


- I 

-1  ^ I ^ ^ 

STT^'juje  1 n»l 


Ar. 


■1  ( r , r , r + A r ) 

pc  p cn  cn  c 


N +N 
c t 

+ I 

n»N 


'n+1  n 


Ar . 


'^pt^'^p’  ’^t.n+l-N^’  "^t  .n+l-N^"^  '^'^c  ^ 

p • c , t (4  5) 


where 


‘ pq 


"n-* 


dr  ' 


+ cos 


d<J) ' 


p . q =•  c , c 


(46  ) 


R 

pq 


1 

dr  ' 


d<l)  ' 


p . q = c , t 


(47  ) 
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The  inner  integrals  in  (46)  and  (47)  can  be  approximately 
analytically  integrated.  They  are  all  of  the  form 


(4  8) 


where  R is  of  the  form 

R = ’^r'  2r  ' b + c (4  9) 

Since  the  range  of  integration  + ^ small  com- 

, - j kR  . 

pared  to  a wavelength,  it  is  appropriate  to  expand  e in 

a Taylor  series  about  some  point  R^  which  is  the  distance 

from  the  observation  point  to  a point  r^  in  the  Interval 

[r  , r ] . Thus , 
n~  n-r 

-jkR  -jk(R-R^) 
e = e e 


-jkR 

"[cos  k (R-R^^  ■ j k(R-R^)] 


-jkR  k”(R-R  ) 

"[1 ^ ^ 


jk(R-R^) 


+ jk  ^ (R-R  ) ^/6  ] 
n 


The  error  in  the  real  and  imaginary  parts  is  less  than 


Ma  X 

r ' fc  [ r - , r ] 
‘ n n+  ‘ 


4 4 

k^(R-R  ) 
n_ 

24 


4 4 

k^(Ar/2)^ 

24 
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where  Ar  is  the  subdomain  size.  For  five  subdomains  per 
wavelength  (Ar/A  = 1/5),  this  results  in  a maximum  error 
of  less  than  1%  in  both  the  real  and  imaginary  parts  of 
the  integrand.  The  resulting  integral  should  indeed  be 
much  more  accurate  than  this.  With  this  approximation, 


I 

/ 


n+  -JkR 


whe  r e 


.■/  ^ 


’-k^(R-R  )^/2 
n 


dr  ' 


k^R^ 

R 

r + b 

1- 

n 

Zn 

n + 

n + 

2 

R + 

r + b 

. 

n- 

n- 

+ k R (r  r^  ) 

n n+  n - 


^r  +b  r+b  „k- 

, 2 r n+  „ n-  „ , c-b  „ . 

-k  [ 2“  ^n+ 2 =^n-  + — 


R +r  +b 
n+  n+ 


R +r  +b 
n-  n- 


and  where 
r 


k(R-R  ) - k^(R-R  )^/6 

^2“  I ' R ' 


I 


(-k  R + k^R^/6)f,n 
n n 


R + r . + b 
n+  n + 


R + r + b 
n-  n- 


+ (k-k  R“/2 ) (r  , - r ) 
n n+  n - 
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k^R  rr  .+  b 


r r . 

n n+ 

2 


r + b . 2 

R ^ R + In 

n+  2 n - 2 


R + r + b 
n+  n+ 


R + r + b 
n-  n- 


, 3 r j.  ~ ^ 

k n+  n- 

~ 6 L 3 


b(r  ? - r ^)  + c(r  - r 
n+  n-  n+ 


.-0 


Combining  these  results,  one  has,  finally, 


I 

/ 


n+  -jkR 

.5— dr'  = 


-jkR„lr  ^2r  2 

1 + jkR^ ^ - j 


3 3 

k R 

n 


( 2 

k^  R 

f 2 i 

R + 
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c-b 

In 
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n+ 

1 2 
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n- 
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T 3 3 

+ jjS _2+ ELZ-. 
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r 9 R 1,2  -I 

+ j-jk  + k R^+  j— 2 + ('^n  + " ’^n-^ 


(5  0) 


Equation  (50)  Is  singular  in  }»'  if  the  observation  point  is 


in  the  interval  [r  ,r  Hence,  the  integrals  in  (-^  6 ) and 

'■  n-  n + 
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(47)  need  to  be  evaluated  by  subtracting  the  singularity 

from  the  integrand  and  adding  its  Integral  to  the  numerically 

determined  integral.  If  r is  in  the  interior  of  the  inter- 
® n 

val,  r <r  <r  , (5  0)  behaves  like  -22,n|'})'|  near  O'  = 0; 

n—  n n + 

if  r =•  r or  r = r . (50)  behaves  like  -£n  I 0 ' | . The  de- 

n n-  n n+ 

tails  of  the  procedure  parallels  that  described  at  the  end 
of  Section  III. 


SECTION  V 


NUMERICAL  RESULTS  AND  CONCLUSIONS 

This  section  describes  numerical  results  obtained 
from  the  computer  code  developed  from  the  theory  described 
in  Section  IV.  The  resulting  code  was  written  to  model  the 
cone  with  or  without  a topcap.  Hence,  results  from  this 
general  code  could  be  checked  against  those  obtained  from 
the  code  based  on  the  methods  of  Section  III  for  the  cone 
without  a topcap.  For  narrow  cone  angles,  the  calculated 
input  impedance  for  unloaded  cones  for  various  frequencies 
was  also  compared  to  the  theory  of  Schelkunoff  [5]  and 
found  to  be  in  very  good  agreement.  Results  from  the  general 
code  for  moderate  cone  angles  were  also  compared  with  those 
computed  by  the  method  of  Appendix  B,  which  includes  the 
effects  of  the  topcap.  These  comparisons  were  made  to  vali- 
date the  consistency  of  the  various  approaches  and  to  com- 
pare with  existing  data.  It  was  also  established  that 
the  input  reactance  at  low  frequencies  could  be  used  to 
check  the  static  capacitance  calculated  in  the  companion 
report  [1]  for  both  the  loaded  and  unloaded  case.  Finally, 
it  was  verified  that  the  computed  results  were  almost 
Independent  of  the  choice  of  the  waist  radius,  a,  of  Section 
III,  provided  a was  chosen  small  enough. 
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All  the  data  in  this  section  pertain  to  a conical 
antenna  with  a vertical  height  of  40  meters  and  a cone 
angle  of  9^“  42.26°.  These  parameters  translate  to  a cone 
slant  height  of  54.05  meters  and  correspond  approximately 
to  the  cone  considered  in  [6].  The  locations  and  values 
of  the  lumped  resistive  loads  used  are  listed  in  Table  1 
and  are  taken  from  [6]. 

Figures  7-10  illustrate  the  current  distribution  on 
the  cone  at  a frequency  of  825  KHz,  approximately  the  first 
resonant  frequency  of  the  unloaded  structure.  Figures  11-14 
illustrate  the  same  results  at  1.375  MHz,  approximately 
half-way  between  first  and  second  resonance  of  the  unloaded 
structure  (see  Figures  15  and  16).  Two  features  of  the  cur- 
rent distributions  are  notable.  First,  the  edge  condition 
[7],  which  requires  chat  the  current  at  the  edge  has  in- 
finite slope,  and  the  continuity  equation  relating  current 
and  charge,  which  requires  that  Che  total  current  approach 
zero  with  zero  slope  at  the  center  of  the  topcap,  combine 
to  limit  the  amount  of  current  Che  topcap  can  support. 
Secondly,  the  loading,  which  Increases  to  a maximum  at  the 
edge,  further  limits  current  flow  on  the  topcap. 

Figures  15-18  illustrate  the  variation  with  frequency 
of  the  Input  Impedance  of  the  conical  structure  for  the  vari- 
ous loading  and  topcap  configurations.  Again,  the  influence 
of  the  topcap  is  found  to  be  negligible.  The  absence  of 
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Table  1.  Positions  and  values  of  loading 
resistors  on  the  cone. 


LENGTH  ALONG 
CONE  GENERATOR 

RESISTANCE 

(METERS) 

(OHMS) 

12.17 

4 .69 

14.34 

7.17 

16.90 

9 .06 

19.  93 

11.74 

23.  54 

15.43 

27.73 

20.82 

32.  59 

30.36 

38.41 

50. 33 

45.31 

114.88 

53.43 

114 . 88 

63.15 

100.00 

72.25 

100.00 

81.35 

100 . 00 
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Figure  12.  Current  on  unloaded  cone  without  topcap, 
L=54.05m,  9q-  42.26°,  V - 1 Volt, 
f » 1.375  MHz . 
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resonances  in  the  loaded  case  can  be  attributed  to  the 
effectiveness  of  the  resistive  loading  in  eliminating 
reflections  from  the  cone  edge  which  would  result  in 
standing  waves  on  the  structure.  However,  one  would  also 

expect  the  bicone  input  impedance  to  approach  a value  of 
56.5  + j 0.0  ohms,  the  input  impedance  of  an  infinite  bicone 
of  the  same  cone  angle,  with  increasing  frequency.  Instead,  the 
impedance  in  Figures  17  and  18  seems  to  be  approaching  a value 
nearer  AO  + j 0.0  ohms.  Computations  at  higher  frequencies 
indicate  that  the  real  part  of  the  impedance  begins  to  increase 
at  about  3.0  MHz  and  at  10  MHz  is  at  about  60  + j 8 ohms, 
slightly  above  the  infinite  bicone  impedance.  Repeating  the 
computations  with  the  cont inuous  loading  function  described 
in  [6]  did  indeed  result  in  an  input  impedance  which  mono- 
tonically  approached  that  of  an  infinite  bicone.  Thus,  it 
appears  that  some  further  adjustment  in  the  valu  of  the  dis- 
crete loading  resistors  might  be  made  in  order  to  more  closely 
approximate  the  impedance  behavior  of  the  infinite  bicone. 

Radiation  patterns  in  the  near-field  region  of  the 
loaded  structure  with  a topcap  (r  • 100  meters  and  the 
frequency  is  550  KHz)  are  shown  in  Figures  19-21.  Figures 
22-2A  give  the  corresponding  patterns  in  the  far  field 

4 

(r  » 10  meters).  For  comparison,  far  field  patterns  tor 
the  unloaded  structure  are  illustrated  in  Figures  25-27. 

Although  resources  did  not  permit  a time-domain  analysis 
of  the  response  of  the  structure,  such  a study,  which  could 

; include  a simple  equivalent  circuit  model  of  the  pulser, 
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Figure  21.  E radiation  pattern  for  loaded  cone 
with  topcap,  L=54.05m,  9 = 42.26°, 
V„-  1 Volt,  f - 550  KHz.°r  - 100m. 
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with  topcap,  L»54.05m,  9^*  42.2 
Vq»  1 Volt,  f - 550  KHz,  r - 10 


54 


would  be  a logical  extension  of  the  present  problem.  To 
be  done  efficiently,  however,  some  improvements  in  the 
present  computer  code  should  be  implemented.  Specifically, 
an  adaptive  integration  procedure  should  be  employed  to 
handle  the  integrations  over  the  conical  current  sub- 
domains, whose  radii  vary  drastically  from  regions  near  the 
feed  to  those  near  the  cone  edge.  The  present  code  uses  a 
fixed  order  quadrature  rule  for  all  segments  on  the 
structure.  Additional  parameter  studies  can  be  carried  out 
using  the  present  code  to  assess  the  effects  of  lumped  v s . 
distributed  loading  and  the  effects  of  various  load  distri- 
butions on  the  performance  of  the  simulator  A more 
ambitious  project  would  more  carefully  model  the  actual 
wire  structure  with  loading. 

One  conclusion  of  this  and  the  companion  study  [l] 
is  that  the  addition  of  a topcap  does  not  significantly 
change  the  electromagnetic  parameters  of  the  structure  - 
at  low  frequencies,  the  static  capacitance  and  effective 
heights  are  almost  unchanged  and  at  the  higher  frequencies, 
the  loading  and  the  sharp  angle  at  the  edge  tend  to  pre- 
vent current  from  flowing  on  the  topcap.  This  observation 
may  have  some  impact  on  the  design  of  future  simulators. 
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APPENDIX  A 


CALCULATION  OF  RADIATED  FIELDS 


Once  numerical  values  for  the  current  distribution 
have  been  determined,  the  fields  radiated  by  the  bicone 
structure  can  readily  be  determined.  Because  of  the 
symmetry  of  the  structure  and  the  excitation,  the  only 
non-zero  components  of  the  electric  and  magnetic  fields 
are  E , E, , and  H,.  These  are  defined  in  terms  of  the 

T 

vector  and  scalar  potentials  as 


E = -i(jjA 
r r 


3r 


■ j WA, 


1 3<P 

r re 


UH, 


-A 
r 0 


3 A. 

c 

3 r 


, 3A 

i : 

r 39 


(A-1) 


where  A « A r + Aq9.  A spherical  coordinate  svstem  centered 
r o 

at  the  bicone  feed  and  with  9 measured  from  the  z-axis  is 
assumed.  Since  the  fields  are  (^-independent,  all  fields 
are  evaluated  in  the  x-z  plane  where  $ = 0.  The  vector 
potential  and  scalar  potential  are  given  by 

I , 


A (r,9) 
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8-^' 
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T ’"l 


+ I IT  (r,e,r)+ 

. n pc  n 

n«  2 
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'N  +1 
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4 ^ ■*■  *^pt  ® ’ ’^N^  + 1 ~ 4 


+ I I4'  (r,9,r  )>;p  = r,9 

^..  . - n p t n I 


n=N  +2 
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(A-2) 


$(r,9)  = ^4—  I pjH'  K 
SiTjuje  n=l  - c 


N I _^,-I  Ar^ 

-I 

n=N  +l  I ^ 

c 


(A-3  ) 


where  che  currents  and  coordinates  are  defined  in  Section  IV. 
The  potential  functions  4'  and  ip  are  defined  as 


4'  (r,9,r  )-Ar 

pq  q q 


e , r~ 

cos  C — ; + cos  ^ 

DO  _+  t 


dO'  (A-4) 


■p  (r,9,r  )=Ar 

q q q 


d(j)',  p=r,9;  q = c,t  (A- 5) 


where 

cos  £ ~ » ± sin  9 cos  $'  sin  9.  + cos  9 cos  9. 
rc  0 u 

cos  ^0^  * - cos  9 cos  i})'  sin  9^  - sin  9 cos  9^ 

cos  C ” * + sin  9 cos  <p’ 
rt 

i “ 

cos  ^0^  * cos  9 cos  (J) 

and  the  radius  vectors  are  all  of  the  form 

+ . m , I ■ 

O “ J 0 ♦ > 

*^0  ^r’  + 2b'r'-*-c' 

q q q q 


with 


» -r  cos  <j) ' sin  9 sin  9^  + r cos  9 cos  9^-3  cot  9^ 

±2  2 2 

c = r ± 2ra  cos  9 esc  9.  + a esc  9„ 
c 0 0 

+ 

= -r  sin9cos  0' 

± 2 - 2 

c » r + 2rL  cos  9 cos  9.  + L 
c 0 

Assuming  a suitable  choice  for  Ar  and  A9 , one  may  approxi- 
mately compute  the  fields  in  (A-1)  by  finite  difference 
approximations ; 


E^(r  ,9) 


, , . ^ 1 ^ , 9)  -<t>  (r  , 9 ) 

-jo)A^(r,9) — 


E-(r,9)  > -ju)A^(r,9)  - ^ ^ ^ ^ 


rA0 


, Ap ( r+Ar ,9 ) -A- (r , 9 ) 

UH.  - -A.(r,9)  + — 2 

'})  r 9 Ar 


A^(r,9+A9)-A^(r,9) 

rA9 


(A-6  ) 
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APPENDIX  B 


AN  ALTERNATE  INTEGRAL  EQUATION  FOR 
A CONE  WITH  TOPCAP 


The  purpose  of  this  appendix  is  to  show  how  a novel 
identity  involving  the  free  space  Green's  function  may  be 
used  to  change  the  integral  equation  into  a form  where  the 
testing  procedure  of  Section  III  is  applicable.  In  ex- 
change for  simplicity  in  the  form  resulting  from  the  testing 
procedure,  however,  one  obtains  extremely  complicated  ker- 
nels in  the  integral  equation.  Furthermore,  the  new  ker- 
nels have  a number  of  singularities  other  than  the  usual 
one  where  match  points  and  field  points  coincide.  These 
complications  make  both  the  analysis  and  the  numerical 
treatment  tedious.  Nevertheless,  numerical  results  have 
been  obtained  for  several  cases  using  the  approach  and  the 
results  are  in  good  agreement  with  data  obtained  by  the 
method  of  Section  III.  For  simplicity,  we  treat  here  only 
the  unloaded  cone. 

As  a prelude  to  the  integral  equation  derivation, 
we  derive  a transformation  of  the  formula  for  electric 
field  components.  Consider  the  x-coraponent  of  electric 
field  given  by 


jfjjue  E^ 


X- (k  +VV* ) A 


k^A  + ;^(V-A) 

X 3 X ' ^ 

^2  . 3 ‘A 

( X + k ) A *f  r — ^ 

3x^  X 3x3y 
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3^A 
2 

3 x3  2 


(B-1) 


where  the  vector  potential  A in  terms  of  current  density 


/-j  K.|  r-r 

I r-r  ' 


The  identity 
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1-  -f4  * 


2_^  2 R 
V +w  [ 


/ 2 2 2 

where  R « u +v  +w  can  be  used  with  u = x-x' , v 
and  w = z-z'  to  rewrite  (B-1)  as 


j'ijue  E 


(x-x’) (y-y') 


r (x-x  ' ) ( 2-Z  ' ) 

^ (y-y ' ) ^+(z-z ' ) ^ 


-jkR 


• X - 


(x-x'  ) 


' \ rn  /■ I 


z-z’  ) 


(y-y’)2+  (z-z’)^ 


-j  kR 


X ^ dV 


as  the  sum 


where 


r-r'  = (r-r’)^x  + (r-r’)^£ 


(r-r’ - X* (r-r’ ) - x-x’ 


is  just  the  component  of  r-r'  along  the  direction  of  x and 
(r-r')^£  = r-r'  - (r-r')^x 

is  just  the  component  of  r-r'  transverse  to  x.  Thus, 

(B-4)  can  be  written  as 


£(?-r ')(?-?') 

t x 


(?-r  ' ) 


2 

t 


R 


dV  ' 


(B-5) 


Since  the  choice  of  the  coordinate  system  is  arbitrary,  we 
may  choose  the  x-axis  parallel  to  some  constant  unit  vector 
a and  write  the  component  of  electric  field  in  the  direction 
of  a to  be 


jcjue 


t(r-r')|.(r-r')a 


(r-?'  ): 


X 


- j kR 

dV 


£ cot 


(a , r-r 


dV 


(B-6) 


where  now  £ denotes  the  direction  of  the  component  of 
(r-r')  transverse  to  a and  s denotes  distance  along  a line 
in  the  direction  of  a.  Note  that  the  integrand  in  (B-6)  is 
singular  not  only  when  R»0,  but  also  when  the  angle  between 
a and  r-r'  becomes  either  0°  or  180°.  Since  in  (B-6)  the 
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only  differential  operator  is  the  harmonic  operator,  then 
along  a line  in  the  direction  of  a,  the  testing  procedure 
of  Section  III  which  uses  piecewise  sinusoids  may  again  be 
used  to  transform  the  harmonic  operator  into  a finite  dif- 
ference operator.  This  is  the  advantage,  gained  at  the 
expense  of  obtaining  a more  complicated  kernel,  of  em- 
ploying the  transformation  (B-3). 

Returning  to  the  cone  problem,  we  choose  the  direction 
of  a to  be  along  the  cone  generator  formed  by  the  inter- 
section of  the  (})"0  plane  and  the  cone  surface,  and  apply 
the  boundary  conditions.  After  some  s t ra igh t f o rwa rd  but 
tedious  vector  projection  operations,  one  arrives  at  the 
integral  equations 
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The 


kernels 


K“  in  (B-9)  are  of  Che  form 

pq 


K~  (r 

pq  p q 


^ C-  (r',^')[A-  ((()’)r'+B-  (Ct.’)l 

D-  (4,-)+-E5_3 E3 a_2i 
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pq 


pq 


The  distance  between  source  points  and  field  points,  ^pq» 
is  of  the  form 


+ y 2 + ! ± 

R~  - "^r  ' + 2b”  r'  + c 

pq  q pq  q pq 


and  b”  and  c”  are  as  defined  in  Section  IV  with  "a"  set 

pq  pq 

+ 

equal  to  zero.  The  term  C~  may  be  expressed  as 

^ pq 


■*•  + + 

C”  = e“  r'  + f” 

pq  pq  q pq 


+ + + + 

The  parameters  A , B , D”  , and  f”  are  defined  in 

pq  pq  pq  pq 

Tables  (B-l)-(B-4). 


Testing  Equations  (B-7)and  (B-8)  with  piecewise 


sinusoidal  testing  functions  as  in  Section  III  results  in 


the  equations 


jwues inkAr 


+ l’cc<'c2>  * -'o 


(B-10) 
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TABLE  B-1 


DEFINITIONS  OF  PARAMETERS  FOR  THE  KERNEL  K 


c c 
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sin*"  6-[sin^(J)'  + cos^9_(cos  <t> ' + 1)^] 
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- 0 
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- e 
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TABLE  B-2 
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TABLE  B-3 


DEFINITIONS  OF  PARAMETERS  FOR  THE  KERNEL  K 


sin^9g  sin^4>'  + cos^Qq 


L cos  0 , 


-e  = sin  9„  cos  ft' 
tc  0 


TABLE  B-4 


DEFINITIONS  OF  PARAMETERS  FOR  THE  KERNEL  K 
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Testing  at  the  cone  edge  with  a piecewise  sinusoidal  testing 
function  which  straddles  both  the  cone  and  the  topcap  and 
which  has  its  peak  value  at  the  cone  edge,  one  obtains 

— r-T l-coskAr['l'  (r  )+4'  (r.,  )] 
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+ ^ 'i'  ( r „ , ) + 4'  ( r „ , ) ]\ 

cc  c,N  -1  ct  c,N  -1  ‘ ] 

c ’ c J 

+ ^ T— T7 I-  COS  kAr  [4'  (r  ,)  + 4'  (r  , ) ] 
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Finally,  testing  on  the  topcap  surface  yields 
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Substitution  of  the  current  expansions,  Eqs . (40)  and  (41)  of 

Section  IV,  into  ( B- 1 0 ) - ( B- 1 3 ) y ie Id s a matrix  equation  for 
the  determination  of  the  unknown  current  coefficients.  Be- 
cause of  the  pulse  expansion  for  the  current,  the  matrix 
elements  involve  integrals  like  (B-9)  but  with  the  current 
in  (B-9)  equal  to  unity  and  the  limits  on  the  radial  inte- 
gration replaced  by  the  limits  of  the  corresponding  current 
subdomain.  According  to  (B-12),  the  term  at  the  edge  also 
requires  the  derivative  of  such  integrals.  In  the  following, 
we  present  a procedure  for  approximately  evaluating  the 
radial  integration,  leaving  the  0'  integration  to  be  done 
numerically.  The  required  integrals  are  all  of  the  form 
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where,  for  convenience,  all  subscripts  and  superscripts  have 

been  suppressed.  Since  the  number  of  subdomains  should  be 

chosen  such  that  k|r^_^-  r^  | is  small,  we  chose  some  point 

r^  in  the  interval  t'^r,+  ’'^n  ^ expand  exp(-jkR)  in  a Taylor 

series  about  the  point  r'«  r ; 

n 
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-jkR  -jkR  -jk(R-R  ) 
e = e e 
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(B-15) 


where  R denotes  R evaluated  at  r'»r  . 
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The  resulting  approxi- 


mate integral  is 
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The  first  integral  may  be  evaluated  by  substituting  the 
definitions  for  R and  C in  terms  of  b,c,e,  and  f; 
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evaluated  at  r'=r  , r ,,  and  r , respectively. 

n n+  n- 

The  evaluation  of  is  facilitated  by  expanding  it 

in  partial  fractions  and  using  the  substitution 
/ 2 

sinhd  = (r'+b)/  c-b  to  obtain 
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/ 2 / 2 
c-b  cosh  9 - e c-b^sinh  9 - f+be 
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ing  CR  2.451.2  and  GR  2.451.4,  one  finds  the  latter  integral 
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Equations  ( B- 1 7 ) - ( B- 1 9 ) complete  the  evaluation  of  (B-16). 


The  derivative  terms  appearing  in  (B-12)  require 
evaluation  of  integrals  of  the  form 


as  appropriate. 


. 


where  the  unprimed  variable  r is  r^  or  r^. 

The  two  edge  terms  also  have  non- in t eg r ab 1 e singularities 
at  the  edge  which  cancel  between  the  various  terms.  To 
handle  this  situation  numerically,  the  singularity  must  be 
explicitly  identified  and  removed  for  numerical  integration. 
Thus  the  same  kind  of  approximate  analytical  integration  of 
(B-20)  as  used  to  evaluate  (B-16)  would  both  eliminate  one 
integration  and  explicitly  identify  the  singular  term.  The 
derivative  can  be  taken  inside  the  integral  if  care  is  taken 
to  identify  the  singular  terms.  Noting  that  9C/3r  = 1, 

3 (R^-C^)/3r  * 0,  and  3A/Dr  = 3B/3r  = 0,  we  have 


/n+  _ n+ 

K...  -J-  [ 


-jkPC  DC 


n- 


n- 


Ar  ’ +B 

2 2 
R -C 


(L  _ ikci  . £1  ] 
U r2  R^j 


e dr 


With  the  approximation  of  (B-15),  the  above  may  be  written 

- j kR 


/ If-' 


= e 


+ I3  + + I?  + I3)  (B-21) 


n- 


The  various  iategrals  appearing  in  (B-21)  are  defined  and 
evaluated  as  follows: 
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¥ 


n+ 


D(er '+f ) 


dr  ' 


• 

R +r 

■* 

De(R  ,-R  )-(be-f)2n 

n+  n- 

n+ 

n+ 

R +r 

+ b 

n- 

n- 

- 

(B-22) 


where  Dw  380.001  and  Dw  380.011  have  been  used. 


= k R 


1 (De-A)r  +f 

"J  r'^  + 2br' 


+ fD-B 


dr  ' 


+ c 


!n± 

2 

= k R 

(De-A)  J.n 

n 

R 

n- 

f D-B-b (De 

-A) 

/ 

c - b 


Can 


1 .L+b 
-1  n+ 


- tan 


/ ,2 

c -b 


, r +b 
- 1 n- 

/ ^ 

c -b 


where  Dw  160.01  and  Dw  160.11  have  been  used.  L’sing 


(D-23) 


Dw  300.003  and  Dw  380.013,  one  obtains 

r 


n+ 


= (1+jkR 
o r 


l+jkR, 


1 

u r 

J ^ 

r 

n- 

r 

f r r 

( [B-Df-b(A-De) ] 

n+  n- 

R R 

1 

1 n+  n-^ 

c-b 


+ [b  (B-Df )-c (A-De) ] 


• n+  n-’' J 


(B-24  ) 
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Again  using  Dw  160.01  and  Dw  160.11,  we  have 


I7  = -jk(l+jkR 


1 


n+ 


Ar  '+B 
R -C 


dr  ' 


• j k ( 1+ j kR 


•'/ 


n+ 


Ar  ' +B 


dr  ' 


(l-e^)r’ ^+2(b-ef )r '+c-f^ 


-jk(l+jkR^)  ^ 


2(l-e") 


iln 


2 2 
R _f-C  ^ 
n+  n + 

2 2 
R -C 
n-  n- 


^ B(l-e‘^)-(b-ef  )A 


tan 


-1 


^(1-e  )r  +b-ef] 
n+ 


- can 


-1 


(1-e  )r  +b-ef 
n- 


J 

(3-25) 


where 


- ’^d-e")  (c-f^)-(b-ef)^ 


The  remaining  integral  is 
r 


Ig  = k 


I 


n+ 


(Ar  *+B) 


-k‘ 


n- 

/ 


R ^2j 

R -C  ‘ 


n+ 


C .(Ar’+B)  , 

2 2 

R(R  -C  ) 
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* -k^ 


r 

/ 


n+ 


C(er ’+f) (Ar ’+B)  , 

2 2 

R(R  -C  ) 


-k 


/ 


n + 


2 f CfAer '^+(Af+Be)r'+Bf , 


2 2 
R(R  -C  ) 


Dividing  R^-C^  “ (l-e^)r'^+  2 ( b- e f ) r ' +c - f ^ into  the 
bracketed  term  in  the  numerator  of  the  integrand,  we  may 
write  the  integral  as 


^8  = 


r £ 

Ae 

. (Af+Be+W) r '+Bf+U~ 

1 ' 

, 2 

.2  ^2 

1-e 

R -C 

dr  ' 


where 


2Ae (b-ef ) 
l-e^ 


Ae(c-f 

1 2 
1-e 

Expanding  the  second  term  in  brackets  in  the  integrand  in 
partial  fractions,  one  can  write  I„  as 

o 

I - I ' + I " 

^8  8 8 


whe  r e 


^8’ 


-k^Ae 
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1-e 


1 

/ 


n+ 


er  '+f 


dr  ' 


-k  Ae 
1 - e “ 


R +r 

. +b 

— 

e(R  ^-R  ) + (f-be)in 

n+  n- 

n + 

n+ 

R +r 

+ b 

n - 

n- 

I 

i 
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(B-26) 


and 


-k‘ 


r 

J - -jL 


dr  ' 


n- 

/ 2 

The  substitution  sinh  9 = (r'+  b)/  c-b  enables  one  to  write 
Ig"  as 


-k 


n+ 


/ 2 

[(Af+Be+W)  c-b  sinh  9 + Bf+U-b ( Af+Be+W) ] 


n- 


/ 2 / ^ 
c-b  cosh9-e  c-b  sinh9-f+be 


/ 2 v''  2 

c-b  cosh  9+e  c-b  sinh  6 + f-be 


d9 


Using  GR  2.451.2  and  GR  2.451.4,  we  obtain  finally 


, 2 Af+Be+W  „ 
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- tan 


'■  '^c-b‘  (R  +C  J + (be-f)R  +b^-c 
n+  a+  n+ 


+ tan 


r / 2 

1 c-b  (R  +C  ) + 
n-  n - 

(be-f)R  +b^-c' 

n- 

+ b)F, 

1 n- 

5 

> 

(B-27) 


where 


[Bf-t-U-b  (Af  + Be-KJ)  ] ( 1- e ) - e ( be- f ) (Af+Be-t-W) 
(l-e")F, 


Equations  (3-22)  through  (B-27)  complete  the  evaluation  of 

the  integral,  (B-21).  Recall  that  the  integral  (B-21)  needs 

to  be  evaluated  only  for  observation  points  at  the  bicone 

edge  (see  Eq . (B-12)).  For  the  source  current  pulse 

associated  with  the  bicone  edge,  there  results  a non- 

integrable  singularity  (with  respect  to  (}) ' integration) 

which  comes  from  the  term  1/R  , in  (B-24).  Each  of  the 

n + 

derivative  terms  in  (B-12)  contains  such  a non-lntegrable 
singularity,  however,  and  they  are  of  opposite  signs  so  as 
to  cancel  each  other.  For  numerical  integration,  of  course, 
the  canceling  singularities  must  be  analytically  subtracted. 

The  integrals  I through  Ig  contain  integrable  singu- 
larities such  as  the  usual  one  where  source  and  field  points 
coincide  (i.e.,  R=0).  In  addition,  however,  there  are  also 

integrable  singularities  introduced  by  the  transformation 
(B-6).  These  arise  from  current  sources  which  lie  along  and 
are  directed  transverse  to  the  line  which  passes  through 
the  observation  point  and  which  is  in  the  direction  of  the 
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electric  field  component  of  interest.  For  example,  in 

computing  the  tangential  field  component  E on  the  cone, 

c 

singularities  arise  from  currents  on  the  topcap  and  its 
i'.iage.  While  the  many  singularities  complicate  the  numerical 
procedure,  they  can,  in  principle,  be  handled.  However, 
the  unwieldiness  of  the  functions  appearing  in  the  integrals 
through  Ig  makes  the  numerical  procedure  rather  ineffi- 
cient and  subject  to  error.  For  example,  the  computer  code 
derived  from  the  formulation  presented  here  seemed  to  yield 
reasonable  results  for  moderate  cone  angles,  but  often 
yielded  erroneous  results  for  the  very  small  cone  angles 
used  to  check  the  program.  The  complexity  of  the  formulation 
became  a considerable  hinderance  in  determining  the  source 
of  these  difficulties.  Consequently,  the  final  calculations 
were  done  using  the  formulation  of  Section  III  which  was 
developed  as  an  extension  of  methods  currently  being  used 
to  treat  flat  plate  surfaces. 
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